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Abstract 

From string theory and the observation of a positive vacuum energy in our uni- 
verse it seems inevitable that there will eventually be a phase transition to an exactly 
supersymmetric (susy) universe. In this phase there will be an effective weakening of 
the Pauli principle due to fermi-bose degeneracy. As a consequence molecular binding 
will be significantly affected. We make some general comments on susy molecules and 
perform a variational principle estimate of ionic binding energies. 



1 Introduction 



Very little is known at present about the behavior of bulk matter in a supersymmetric 
background. Since supersymmetry is a prominent feature of string vacua a prime question 
is whether or not the evolution of intelligent supersymmetric life forms would be possible 
in such a background. If not, since we could hardly expect to be making observations in a 
universe which did not support life, observer bias would offer some type of explanation for 
why our universe has, or seems to have, a broken supersymmetry (susy). In this sense a 
fuller understanding why our universe is as it is might require an exploration of alternative 
universes. 

A supersymmetric universe, at least in flat space, would be expected to have a lower 
vacuum energy than the observed positive vacuum energy of our current universe. There- 
fore, from the string theory point of view, the current universe is intrinsically unstable [1] 
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to a future exactly supersymmetric universe with, possibly, zero vacuum energy. A first 
phenomenological glance at the properties of such a universe has been already undertaken 
[2]. As discussed there the transition will begin with the quantum nucleation of a bubble 
of true (susy) vacuum of radius sufficiently large to allow bubble growth. It is likely that 
this critical size is of galactic dimensions [3] . As the scale of the universe expands, the prob- 
ability per unit time for the nucleation of a critically sized bubble will grow without limit 
until the transition occurs. If the current acceleration continues until all matter is trapped 
in white dwarfs, neutron stars, and black holes, the eventual phase transition will lead to a 
rapid series of gamma ray bursts [I] followed by collapse of all matter into susy black holes 
[5]. On the other hand, if the transition occurs soon enough that earth-like planets are still 
orbiting stars, susy planets will be formed and it is possible, depending on susy molecular 
physics, that susy life forms will evolve [2]. The modest goal of the current paper, however, 
is solely a study of the properties of ionically bound molecules assuming susy atoms come 
into existence in the form of small, heavy nuclei surrounded by an electron/selectron cloud. 

In our broken susy universe, the structure of matter is dominated by the Pauli principle 
which forces bound fermions into high energy levels. After a transition to the exact susy 
universe, these excited fermions will convert in pairs to their degenerate scalar partners 
through gluino or photino exchange: 



The scalars, /, not being constrained by the Pauli principle, will fall into the ground 
state wave function via, for example, photon emission. Thus, in exact susy, there will 
be no more than two fermions of each type in each bound state. The remaining baryon 
and lepton number will be filled out by sbaryons and sleptons with all stable nuclei and 
atoms containing particles in s-waves only. Where possible without confusion, we will refer 
generically to fermionic and scalar electrons as electrons and to fermionic and scalar nucleons 
as nucleons. In exact susy all particles and all bound states are members of degenerate susy 
multiplets. 

Our variational approach allows an approximate but analytic solution for the ion energies, 
mean radii, and ionic binding energies. For the sake of numerical comparison with our world, 
we evaluate our formulae in the case that the masses of the degenerate susy particles will be 
the masses of the standard model particles in the broken susy world. This is approximately 
what would be expected if the ground state of the universe corresponded to one of the higgs 
extensions of the minimal supersymmetric standard model (MSSM) with all susy breaking 
parameters set to zercfl In these models the coupling constants would also be expected to 
be approximately as they are in the broken susy world with small deviations given by the 
running of the couplings over a slightly larger domain of exact susy. Other models for susy 
breaking in our world, such as gauge mediated models, could lead to similar results as the 
susy breaking is turned off. In the model of radiative breaking of electroweak symmetry, 
the symmetry breakdown of the standard model is due to a squared Higgs mass running to 
negative values. This could continue to be true if exact supersymmetry extends into the low 
energy domain 
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One could also entertain alternative possibilities for the ground state of the universe such 
as the supersymmetric AdS/CFT model or one of the five flat space superstring theories. 
These might qualitatively share some of the features we investigate here. We give analytic 
results for the energies and mean radii, so the consequences of other assumptions about the 
degenerate susy multiplet masses and charges can be readily evaluated. 

In section II we approximate the multi-electron wave function as products of hydrogenic 
wave functions. This wave function is improved by the variational method to take into 
account screening in section III. In section IV we calculate the ionic bonding of some diatomic 
molecules that are prominent in our broken susy world and in section V we treat ionic binding 
of the susy water molecule. We conclude with a summary and some discussion in section VI. 



2 Susy ions 

Working in natural units, U — c — 1, the Hamiltonian for the system of nuclear charge Z 
and N electrons is given by: 

i N N i N -i 

Z '" L i = l i=l i<j I 1 ! L jl 

with e 2 being the fine structure constant, a ~ 1/137, and m being the electron mass. 

From the variational principle, an upper bound on the ground state energy is given by 
the expectation value of the Hamiltonian for any trial wavefunction. 

E gs = (V>|i*W • (2.2) 
The ground state wavefunction for the hydrogen atom is given by: 
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We might take the trial wavefunction for the general ion to be the product of this ground 
state wavefunction for each electron. 

i >= ^ me ^y / \-™^ . (2 .3) 

This is an eigenfunction of the first two terms in eq. 12.11 with eigenvalue 

(H) Hydrogenic = -RooZ 2 N (2.4) 

where Roo is the Rydberg constant 

= me 4 /2 = U.GeV. (2.5) 



This wave function is exact for the hydrogenic atoms (N — 1, Z arbitrary) ignoring nuclear 
effects and relativistic corrections. For N > 2 fermionic electrons, this would not be expected 
to be a good approximation since only two electrons can be in this Is wave function. However, 
in the exact susy world, all of the charged leptons can be in the ground state wave function 
with as many of them as necessary becoming bosons via the pair conversion process of eq. ll.U 
Note that the conversion process occurs generically in susy and does not require R parity 
violation. 

Now, we need to calculate the expectation value of the last term, corresponding to the 
repulsion between the electrons 

N , 

V ee = e 2 J2l f . (2.6) 

We calculate (V ee ) directly by integration. 

(V ee ) = e 2 I ^ I V , - |Wr 3 ...dr N . (2.7) 
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The number of terms in the sum is equal to I I . By symmetry, this integral can be 
rewritten and simplified. 

To evaluate this integral, we set the coordinates relative to r 1; so that the angle between 
i"i and Y2 is simply 82. By the law of cosines, 



|ri — r 2 | = yrf + r\ — 2r\r 2 cos 6*2 • 
We start by evaluating with respect to r2, defining a new integral: 



h = f e -2(me*z)r 2 \ r\^mQ 2 dr 2 dB 2 d^ 2 . (2.9) 

J . I r 2 _|_ r .2 _ 2nr 2 cos ^2 



The 02 integral is just 2ir. The 9 2 integral is a bit more involved: 
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This is 2/ri for r 2 < n and 2/r2 for r 2 > r%. 



{n + r 2 - |n -r 2 |) . (2.10) 



Inserting these back into I 2 , we get the following: 

/ 2 =47r (i P e-' ! ^ z y'4dr 2 + e- 2 ' m ' 2z ' r 'r 2 dr 2 \ (2.11) 

=(S^( 1 -( 1+ ™ , ^)^ ai ) ■ (2 ' 12) 

Returning to (V ee ), we now have: 

{V ee )=e 2 (^j 4(me 2 Z) 3 J™ e- 2me2Zri ri (l - (l + me 2 Zr l ) e ~ 2me2Zr ^ dr x (2.13) 
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Including the hydrogen-like term from eq. 12.41 we have 

(H) = -R 0O ZN(z-^(N-l)^ . (2.14) 

By the variational principle this gives an upper limit to the ground state energy or a 
lower limit to the ground state binding energy. 

UN = Z, 

( H ) = -^(3Z 3 + 5Z 2 ). (2.15) 



3 A Better Approximation 

We can improve this upper bound anticipating that electron shielding of the nuclear charge 
will reduce the coefficient in the wave function exponential. With Z s a dimensionfull param- 
eter to be determined by energy minimization we take the normalized trial wavefunction for 

the N electrons to be: 

N ( [73 \ 
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Proceeding as in section 2, we can write the expectation values 



I)=4.f (f)e— r* = Z. , (3,7) 



(3.18) 



and 

Wi-rjl/ 8 

Thus 

{H} = N?i-Ze>Z,N + N{N - 1) e^Z, 
w 2m s 2 8 

This is a minimum at 

Z s = me 2 (z-^{N-l)} . (3.21) 

Since Z s must be positive in order to have a normalizeable wave function, this equation 
is only meaningful for 
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N<— Z+l . (3.22) 
Putting this back into the expression for (H) gives 

(H) = E(N, Z) = -R^N [z - ^(N - 1))" . (3.23) 

Comparing with eq. 12.141 we see that the ground state energy estimate is reduced by 
the improved wave function providing a decreased upper limit to the actual energy and an 
increased lower limit, —E(N, Z), to the actual binding energy. According to the variational 
principle, further improvements in the trial wave function will yield a binding energy ap- 
proximation approaching the exact solution from below. This variational result was known 
shortly after Schroedinger's original result [6] and improves the accuracy of the theoretical 
energy of the helium-like atoms (N = 2, Z arbitrary) by a factor of three above the result 
based on treating the electron repulsion perturbatively [7]. 

Beyond N = 2 the result is not a useful approximation in the broken susy world due to 
the Pauli principle which forbids more than two electrons from occupying the Is level. With 
the circumvention of the Pauli principle in the exact susy case, we can expect comparable 
accuracy for susy atoms of arbitrary iV as for standard helium. 

The average distance between an electron and the nucleus is 

r(N, Z) = 4v3 [°° r 3 dre- 2Z ° r = = -^-(Z - -(N 
ix Jo 2Z S 16 

except that for iV = we must put r(N, Z) = 0. 

We see that, at fixed Z, this mean distance is increased over the hydrogenic (N = 1) 
case due to the electron repulsion. The size of the neutral atoms (N = Z) decreases rapidly 
for increasing Z unlike the case in the broken susy world as discussed below. However, at 
least for moderate Z, the mean radius of the electron cloud remains large compared to the 
nuclear size. 



(3.19) 
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4 Ionic Binding 



Ionic binding occurs when it is energetically advantageous to transfer one or more electrons 
from one atom to another after which the positive and negative ions strongly attract each 
other. Ionic binding, when appreciable, tends to dominate over covalent binding which 
corresponds to two atoms sharing electrons equally. In the broken susy world, ionic binding 
is most prominent in the alkali-halide compounds where one electron can be taken from an 
outer shell of one atom and used to form a closed shell on the other atom. 

The change in energy after transferring / electrons/selectrons from one ion of atomic 
number Z\ to another of atomic number Z 2 and taking into account the subsequent Coulomb 
attraction is 

AE = —E(N\, Z x ) - E(N 2 , Z 2 ) + E{N 2 + I, Z 2 ) + E{N X - I, Z x ) 

e 2 l 2 

(4.25) 



r(N 1 -l,Z 1 )+r(N 2 + l,Z 2 ) ' 

The negative of AE is the molecular binding energy. This expression assumes that the 
perturbation of the atomic wave functions away from spherical symmetry can be neglected 
and that the interatomic distance is rougly the sum of the radii of the two ions. These 
approximations are reasonable as a first attempt to estimate molecular properties in an 
exact susy background. 

Assuming only single electron transfer, The ionic molecular binding energy for two neutral 
atoms of atomic numbers Z\ and Z 2 is 



B = +E(Z U Z{) + E(Z 2 , Z 2 ) - E{Z 2 + 1, Z 2 ) - E{Z X -l,Z 
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+ r(Z 2 + l,Z a )+r(Z 1 -l,Z 1 ) • (426) 



This can be written 

B = F(Z 2 ) - I{Z X ) + e 2 /(n + r 2 ) (4.27) 

where the electron affinity of the atom of atomic number Z 2l namely the energy released 
when an electron is added to the atom , is 

F(Z 2 )=E(Z 2 ,Z 2 )-E(Z 2 + 1,Z 2 ) (4.28) 

and the ionization energy / of atom 1, the energy required to remove an electron, is 

I(Z 1 ) = -E(Z 1 ,Z 1 ) + E(Z 1 -1,Z 1 ) . (4.29) 

Which atom donates an electron and which receives is determined by maximizing the binding 
energy. In general, the atom of lower Z is the donor atom. In eq. 14 .2 71 we have estimated the 
equilibrium separation of the two ions as the sum of the mean radii since the overlapping 
electron clouds will inhibit closer approach. Using the sum of the rms radii does not lead to 
significantly different results. 



In fact, however, we find that, when the binding is positive, the maximum binding comes 
when the atom of lesser Z = Z\ donates all its electrons to the other atom of greater atomic 
number Z 2 . This does not result in a violation of eq. 13.221 

The binding energy is then 

B = E(Z 1 , Z x ) + E(Z 2 , Z 2 ) - E(0, Z 1 ) - E(Z 2 + Z u Z 2 ) + °fl (4.30) 

r(Zi + Z 2 , Z 2 ) 

The mean radius of the totally stripped ion, r(0, Z), is taken to be zero. The molecular 
binding energy for hydrogen with an element of atomic number Z can, therefore, be expressed 
as 

B = E{Z,Z)-E(Z + l,Z)-13.6eV + a/r(Z+l,Z) . (4.31) 

This is the electron affinity of the element of atomic number Z minus the hydrogen ionization 
energy plus the negative of the Coulomb energy of attraction of the two ions. 
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Table 1: single electron ionization energies of alkali atoms in eV in the exact susy case using 
the optimum wave function from eq. 13.161 In the fourth column we give the experimental 
ionization energies in our broken susy world. The fifth column gives the predicted mean 
radius of the electron/selectron cloud in Angstroms. 

In table [T] we tabulate the ionization energies of the alkali atoms and the mean atomic 
radii in the susy case to compare with the experimental values in the broken susy case. It is 
seen that, for large Z in the susy case, the ionization energy is significantly greater than the 
experimental values in our broken susy world as would be expected since all the electrons 
and selectrons are in the Is ground state. In the current universe, all atoms have a radius 
of order 1 A° since, for greater Z, the increased attraction to the nucleus is balanced by 
electron repulsion and by the need to put electrons into higher orbitals which have greater 
spatial extent. In exact susy on the contrary, as can be seen from table Q], heavier atoms 
have significantly smaller radii. In table [2, we give the electron affinities for the halides, 
namely the energy released if an extra electron is dropped onto the neutral atom making a 
negatively charged ion. Table [3] gives the estimated binding energies for the alkali halides in 
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Table 2: electron affinities in eV and mean radii in Angstroms of the halide atoms in the 
exact susy case using the trial wave function from eq. 13.161 
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Table 3: Binding energy estimates in eV from eq. l4.30l for the alkali halides in the exact susy 
case using the trial wave function from eq. 13.161 Negative entries correspond to combinations 
with no ionic binding in the susy world although they are bound in the broken susy world. 
Experimental binding energies in the broken susy case, given in parentheses, are taken from 



the susy case following eq. 14.301 The negative entries correspond to molecules that are not 
bound in the susy case although they do form diatomic molecules in the broken susy case. 
Covalent binding of alkali-halides is negligible in the broken susy world and presumably also 
negligible in the exact susy world at least for those molecules that have significant positive 
ionic binding energy. 

Where the susy binding energies are positive they are much greater than in broken susy 
and the interatomic radii are much smaller. The smaller radii of the susy atoms suggests that 
chemical reaction rates will be much slower in the susy world, while the greater molecular 
binding energies suggests that the diatomic molecules which do bind will be much sturdier 
than in the broken susy case since the dissociation energies reach well into the gamma ray 
region. However, these results are moot if, as suggested by the semi-empirical mass formula 
[2], nuclei above oxygen are rare or nonexistent in a susy background. If there were non- 
negligible surviving heavy elements, due to the very large binding energies in table [3j they 
would rapidly steal atoms from lighter molecules. 

It is also interesting to look at the low-lying diatomic molecules whose susy binding 
energies in the current approximation are given in table |U Unlike the case in the broken 
susy world, in exact susy there is significant binding of helium to other low-lying elements. 
Molecules with small or negative ionic binding energy could be bound through covalent 
binding which is not treated in this article. We find no ionic binding in diatomic molecules 
whose atoms differ by zero or one unit of atomic number except for a few of the lighter 
molecules such as helium hydride which is only marginally bound in broken susy. 

In the case of homonuclear diatomic molecules we can write, following eq. 14.261 the ionic 
binding energy due to the transfer of I electrons from one atom of atomic number Z to the 
other as 



Since this is negative for arbitrary I, the homonuclear diatomic molecules have no ionic 
bonding for the approximate wave functions we are using. Although no ionic bonding is 
found for these potential molecules, the negative binding energies are small compared to 
those of the alkali-halides. This suggests that with a better approximation of the wave 
functions, ionic binding may occur and, additionally, one might expect covalent binding to 
be effective for these molecules. 



5 Supersymmetric water 

As a final example we consider the ionic binding of the susy water molecule. In the current 
(broken susy) universe, water binds due to a transfer of an electron from each of the two 
hydrogen atoms to 2p orbitals of the oxygen atom resulting in a closed shell of oxygen. Since 



B = 2E(Z, Z) - E(Z -l,Z)- E(Z + l,Z) + 
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Table 4: Binding energy estimates in eV from eq. l4.30l for the low lying diatomic molecules in 
the exact susy case using the trial wave function from eq. 13.161 Negative entries correspond 
to unbound combinations. Experimental binding energies [8] in the broken susy world are 
given in parentheses for the molecule indicated by the line above. 



the p orbitals are in one to one correspondence with the Cartesian coordinates, the hydrogen 
bonds, in zeroth approximation, are at right angles to each other. Coulomb repulsion of the 
protons then leads to the slightly greater observed bond angle of 105°. 

In the susy case, all the electrons can be put into the spherically symmetric Is wave 
function of oxygen, converting, if necessary, into scalar electrons via eq. ll.ll Proton repulsion 
will then naturally lead to a linear triatomic molecule with the oxygen ion in the center. Since 
this arrangement has no molecular dipole moment, the properties of water as a solvent might 
be significantly changed in susy matter. 

In this linear configuration, the binding energy of a neutral atom of atomic number Z 
with two hydrogen atoms would be 



128 V ' J y ' 

For H2O this corresponds to a total dissociation energy of 98 eV. Atoms lighter than beryl- 
lium will not have a positive ionic binding energy to two hydrogen atoms. 

6 Conclusions 

The variational approach we have taken in this paper allows for an approximate but analytic 
expression for the atomic energies, mean radii, and ionic molecular binding energies as a 
function of the electron charge and mass. We have tabulated numerical results for the case 
where the degenerate electron/selectron mass is equal to the electron mass in the broken 
susy world but, obviously, other masses can be used if that becomes theoretically motivated. 
Our core assumptions are that, in the susy background as in our world, the electron mass 
is non-zero but much less than the nucleon mass, the nucleus is small compared to atomic 
scales, and atomic structure is governed by a U(l) electromagnetism. The results of the 
paper are also independent of any assumptions about atomic abundances as long as these 
are non-zero. Of course, the degenerate susy masses and the abundances would be important 
parameters in discussing the transition to exact susy and the subsequent properties of the 
susy world. 

Our approximate atomic wave function corresponds to putting all electrons into a screened 
Is wave function. This approximation is of limited usefulness in the broken susy world for 
atoms above helium due to the Pauli principle but could be a good approximation in the 
exact susy world. We have neglected mixing between atoms with the same Z but different 
numbers of scalar electrons. We have computed ionic binding energies of various atoms as- 
suming that atoms remaining spherically symmetric is a reasonable zeroth approximation as 
in the broken susy world. We have also assumed that the bond lengths are approximately 
the sum of the mean atomic radii. This assumption is found to be reasonable in broken susy. 
We have treated the atomic ions quantum mechanically but the molecular bound states 



B = E(Z, Z) - E(Z + 2,Z)- 2^ + 
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are treated classically. A similar semi-classical calculation gives reasonable results for the 
alkali-halides in the broken susy world. 

As expected, we find that the atomic sizes are quite small in general compared to the 
corresponding atomic sizes in the broken susy world. This tends to lead to large ionic binding 
energies since the transfer of an electron from one atom to another leads to a strong Coulomb 
attraction. Such molecules could be expected to be much sturdier in a susy world. On the 
other hand some of the heavier alkali halides that are bound in the broken susy world become 
unbound in exact susy. This is due to the higher ionization energy of the heavy alkali atoms 
in exact susy where all the electrons are in the Is wave function. 

It is interesting to inquire whether the diatomic molecules that are unbound in the susy 
world play any essential roles in the broken susy world. This might lead to an anthropic 
explanation of why our world is not exactly supersymmetric. In fact however, from an 
environmental point of view one readily comes on reasons to welcome the absence of a few 
of the molecules that are unbound in the susy universe but have severe deleterious effects 
in our broken susy world. For example, among the low-lying diatomic molecules, sodium 
fluoride, NaF, potassium chloride, KC1, and potassium fluoride, KF, are known to have life- 
threatening effects in the broken susy world. Of these, NaF and KC1, are seen in table [3] to 
be unbound in the exact susy world. 

The susy water molecule is predicted to be linear rather than forming the familiar 105° 
bond angle seen in the broken susy world. The resulting absence of a dipole moment could 
be significant with respect to viability since, in the broken susy world, the anomalously 
large dipole moment of the water molecule is important in making water a universal solvent 
facilitating chemical reactions. 

Clearly, we have barely begun to explore the properties of bulk susy matter and we are 
far from being able to make any statements about the possibility of life in the susy world. 
In the more proximate future it should be possible to explore the effects of mixing between 
states of differing selectron numbers and to estimate covalent binding energies as well as to 
seek an improved variational wave function. 
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